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The Cygnus Loop is a nearby supernova remnant (SNR) observed across the electromagnetic spectrum. With 
the analysis of 6 years of Fermi/LAT data we find that, what previous studies had considered a single source, 
consists of an extended source plus a point-like source south-east of the SNR. The extended gamma-ray emission 
is well correlated with the thermal X-ray emission of the SNR, and the energy spectrum displays a pronounced 
maximum at ~ 0.6 GeV. However, in a region where the radio emission is strongly and distinctly polarized, the 
gamma-ray spectrum shows no sign of a break. Therefore, the spatially resolved gamma-ray emission permits 
the study of different interaction conditions of the SNR and the surrounding medium. 


1. INTRODUCTION 


The Cygnus Loop is the remnant of a core-collapse 
supe rnova explosion that occurred about 14000 year 
ago Leven son et ah 1998 at a distance of 540lgo° pc 
Blair et ah 2005. The Cygnus Loop is among the 


closest supernova remnants (SNRs) to Earth, which 
implies that it could act as a local accelerator. Due 
to its proximity, the Cygnus Loop is seen on the sky 
with an angular size of about 3 degrees. In general, 
the blast wave of the SNR is not breaking out of a 
dense cloud, but running into a wall of atomic gas re¬ 
lated to the cavity in which the supernova occurred. 
The wall slows down the shock, which becomes bright 
in optical emission lines. The reflected shock prop¬ 
agates through the hot interior, which enhances the 
X-ray em issio n in corr elatio n with the_ o ptica l emis¬ 
sion Graham et al.|199~5| Levenson et al.|1996 . How¬ 
ever, some portions of the shock proceed unimpeded 
through low-density inter-cloud medium. 

X-ray emission from reflection-shocked gas is partic¬ 
ularly bright in the east. In contrast, the south of the 
SNR (the so-called breakout ) is very dim in X-rays. 
This is often regarded as caused by the expansion of 
the blast wave into a low-density medium. However, 


Uyamker et al. 2002 found that the polarization of 


the 2695 MHz emission was much higher there with 
respect to the north of the shell. A possible interpre¬ 
tation of this feature is that a second SNR is present 
in that region, and interacts with the Cygnus Loop. 

No compact object is firmly associated with the col¬ 
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lapsed progenitor of the Cygnus Loop. A few candi¬ 
dates lie within the breakout, where the ASCA survey 
revealed a point-like sou rce, b ut it is n ot firmly estab¬ 
lished as a neutron star Miyata et al. 2001 . There is 
yet another compact object with a candidate pulsar 
wind nebula nearby, r evealed by Suzaku a nd XMM- 

No pul¬ 


Newton observations Katsuda et al 


2012 


sations have been detected from any of these objects. 
In addition, a very high transverse proper motion of 
~ 1300 km s -1 is needed if it is assumed that one of 
these candidate neutron stars departed from the ge¬ 
ometric center of the Cygnus Loop some 14000 years 
ago. Such a supersonic movement would produce a 
cometary shape in the X-ray emission that has not 
been observed so far. However, this could be explained 
if the neutron star was related to the second SNR sug¬ 


gested by Uyamker et al. 2002 


The detection of GeV gamma-ray emission from the 


Cygnus Loop was published in Katagiri et al. 2011 
who analyzed two years of Fermi/LAT data com¬ 
prised between August 2008 and August 2010. In this 
analysis, the shape of the Cygnus Loop was modeled 
as a ring, somewhat more extended than the shell 
seen at other wavelengths. The spectrum is curved 
(modeled as a log-parabola), and the fit to a one-zone 
hadronic model returns plausible values for the pa¬ 
rameters. 


In this work we analyze six years of Fermi/ LAT 
data using the latest software. The factor 3 increase in 
statistics with respect to the previous study provides 
unprecedented sensitivity to study both spatial and 
spectral features of the gamma-ray emission from the 
Cygnus Loop. 
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2. DATA ANALYSIS 

We analyzed Fermi/ LAT Pass 7 Reprocessed data 
corresponding to the period between August 4 t/l 2008 
(start of science operations) and September 7 th 2014. 
We defined the ROI as a circle of 10° radius centered 
at the position (RA, DEC) = (20 /l 58 m ll s , 29 0 23'56"), 
J2000, which is 2° displaced towards negative Galac¬ 
tic latitudes with respect to the catalog position of 
the Cygnus Loop. This is done in order to be less af¬ 
fected by the diffuse emission from the Galactic plane. 
Data were processed with the version v9r32p5 of the 
ScienceTools. We selected class 3 events in the en¬ 
ergy range between 58.5 MeV and 300 GeV, with the 
recommended quality cuts (including the requirement 
for the spacecraft to be in normal operation mode, 
LAT_CONFIG=l, data to be flagged as good quality, 
DATA_QUAL=1, and a cut on the rocking angle of the 
spacecraft, ABS(ROCK_ANGLE)< 52°). In addition, 
we applied a zenith angle cut of 100° in order to pre¬ 
vent event contamination from the Earth limb. Data 
were binned in sky coordinates with the gtbin tool, 
using square bins of 0.125° side. This tool produces 
a counts map , with the number of events recorded by 
the detector. 

We performed a binned likelihood analysis with a 
model containing the standard Galactic and extra- 
galactic diffuse emission models provided in the Sci- 
enceTools , plus the sources in the 2FGL catalog ly¬ 
ing up to 15° away of the ROI center. We call the 
model with the point-like sources plus the Galac¬ 
tic and extra-galactic backgrounds the null hypothe¬ 
sis, which has a maximum likelihood Cq. Then, we 
generate alternative models by adding spatial tem¬ 
plates and by changing the functions describing spec¬ 
tral shape. By varying the parameters of each mod¬ 
els, we compute the corresponding maximum likeli¬ 
hood Cmodei • We choose the best representation of the 
Cygnus Loop as the model which obtains the highest 
value of the likelihood ratio LR=2 log (£o/Anode/)• 

For any of the tested models, we can use the tool 
gtmodel to produce an expected counts map given the 
exposure associated to the data set. For visualization 
purposes, we produce what we call the S/N map by 
subtracting the expected counts map from the actual 
counts maps, and then dividing by the square root of 
the expected counts map. 

Complementary to the Fermi/ LAT data analysis, 
we have re-analyzed the 11 cm radio emission observed 
by_the 100 m Effelsberg telescope. |Uyamker et al 


2002 proposed the Cygnus Loop be divided in two 


regions for the two-SNR interpretation. Based on our 
re-analysis, which considers the presence of extended 
Stokes I radio emission in addition to the distinct and 
intensely polarized radio emission, we have re-defined 
the regions to be more equal in size and both hav¬ 
ing a circular shape. We consider that the southwest 
(SW) feature is the circular region of 1.07° radius cen¬ 


tered at (RA, DEC) = (20 /l 49 m , 29°47') as shown in 
Figure [l] 


3. RESULTS 
3.1. Morphology 

The S/N maps at different energy ranges, produced 
with the null hypothesis are shown in Figure [l] While 
an extended source is clearly seen at energies below 
10 GeV, in the last panel only residual, localized emis¬ 
sion is present south of the SNR. 

As a first s tep w e reproduced the analysis from 
Katagiri et al. 2011| . For this purpose we modeled the 
Cygnus Loop as a ring of 1.6/0.7 degree outer/inner 
diameter, centered at (RA, DEC) = (20 /l 51 m , 30°50'). 
The ring is divided in four quadrants, and the spec¬ 
tral parameters of all of them are varied simultane¬ 
ously. We note that the hard spot remaining above 
10 GeV is included in the southeast (SE) quadrant 
of the ring used in Katagiri et al. 12011]. By sub- 
stituting the SE quadrant by a point-like source, we 
find that LR improves by 124. The position of this 
point-like source optimized by the tool gtfindsrc is 
(RA, DEC) = (20 /l 53 m 55 s , 29°24 / 45 // ) with an uncer¬ 
tainty of 0.02°. We call this source J2053.9+2924. 
Its position is coincident with the X-ray and ra¬ 
dio source 2E 2051.7+2911, which is likely an AGN 
[Brinkmann et ah||1997 . Therefore, we consider that 
J2053.9+2924 is a source in the background of the 
Cygnus Loop, and should not intervene in the model¬ 
ing of the diffuse emission 1 . 

Having included J2053.9+2924 in the list of point- 
like sources, we maximize the likelihood of a template 
generated fro m the X- ray co unts map observed by 
ROSAT Aschenbach and Leahy [1999 , re-binned to 
match the pixel size of maps of the present analysis. 
The likelihood ratios for the spatial models mentioned 
above are shown in Table [H It is clear that the thermal 
X-ray emission correlates very well with the observed 
gamma-ray emission, and requires less degrees of free¬ 
dom than the ring to describe it. 

The X-ray emission is very faint in the region of 
highly polarized radio emission, but we divide the spa¬ 
tial template in order to study this particular region. 
The templates for the main (NE) emission is cropped 
to avoid having pixels accounted for twice in the over¬ 
lapping region (Figure [2]). We verify that both regions 


1 The point-like source found in this analysis is called 
3FGL J2053.9+2922 in the recently publis hed Thir d 
Fermi/LA T Source Catalog The Fermi-LAT Collabora-| 
|tion|2015] . We note that the source overlaps with the template 
for the Cygnu s Loop, which in 3FGL is still modeled as the 
ring defined in Katagiri et al. [2011 . 
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Figure 1: S/N maps of the null hypothesis in different energy ra nges. From l eft to right, the panels include: the two 
regions discussed in s ection 3.2 the ring used for m odeling in|Katagiri et alT| |2Q1 1 ; the radio intensity contours from 


Uyamker et al.| [2002 ; and the X-ray contours from Aschenbach and Leahy | [1999 


Table I Likelihood ratio (LR) of the tested templates, 
with the number of degrees of freedom (d.o.f) added to 
the null hypo thesis after selecting the best spectral model 
(Section 3.2). Three of the additional d.o.f. always 
correspond to the spectral parameters of J2053.9+2924, 


except for (1), where original template by Katagiri et al. 
[2011 is tested. 


Model 

LR 

d.o.f 

0) Null hypothesis 

0 

0 

1) Ring 

Divided in four quadrants 

2069 

12 

2) 3/4 Ring 

SE quadrant substituted by J2053.9+2924 

2193 

12 

3) ROSAT template 

2204 

6 

4) ROSAT template 

Divided in NE and SW 

2238 

8 


contribute significantly to the overall emission. Next, 
we proceed to test different the spectral models. 


3.2. Spectrum 

The energy spectrum of the source is shown in Fig¬ 
ure [3j The global emission, as well as the emission of 
the NE region are well described by a log-parabolic 
shape. To infer the spectral shape of the SW re¬ 
gion, we test a power-law shape, a log-parabolic shape, 
and a power law with exponential cutoff. Because of 
the proximity of J2053+2923 to this region, we also 
test all three possible models for the point-like source. 
Then, we evaluate the likelihood ratio of each combi¬ 


Table II Spectral parameters of the Cygnus Loop and its 
NE and SW regions under the assumption of a 
log-parabolic spectral shape, 

dN/dE = iVo(£/lGeV)“ a “ /31os(£:/1GeV) . The photon 
flux at energies above 58.5 MeV is shown in the last 
column. 


Region 

a 

p 

Flux 

10 8 ph cm -2 s -1 

Cygnus Loop 

2.26 + 0.03 

0.25 + 0.03 

13.5 + 0.9 

NE 

2.24 + 0.04 

0.32 + 0.04 

9.0 ±0.9 

SW 

2.27 + 0.06 

0 

7.2 + 1.1 


nation with respect to the initial assumption of both 
components having power-law shaped spectrum. We 
observe that models where J2053+2923 has an addi¬ 
tional degree of freedom in the spectrum have a like¬ 
lihood ratio with respect to the power-law/power-law 
hypothesis of ~ 10, whereas models where the addi¬ 
tional degree of freedom is added to the SW region 
only improve the likelihood ratio by ~ 1. Therefore, 
we conclude that there is ~ 3 a evidence that the spec¬ 
trum of J2053+2923 is curved, whereas the spectrum 
of the SW region is compatible with being a simple 
power law. The best fit spectral parameters of the 
global emission and the studied regions are shown in 

Table HI 

3.3. The point-like source J2053+2923 

The point-like source South of the Cygnus Loop is 
detected with high significance (TS=214). As men¬ 
tioned in Section |3.2[ the spectrum is described ei- 
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Figure 2: Sketch of the templates used in the Fermi /LAT analysis. From left to right: the total emission, the NE 
region, a nd the SW region. Color code matches that from Figure [3] Green contours are the same X-ray contours as in 
Figure [l] [Asche nbach and Leahy 1999 . 



Figure 3: Energy spectrum of the Cygnus Loop (blue). 
The emission from the NE region is shown in green, while 
the SW region is represented in red. The statistical 
uncertainty range of the best model is shown as a shaded 
area. Spectral points include statistical uncertainties 
(solid bars) and systematic uncertainties (shaded bars). 


ther by a power law with exponential cutoff or by a 
log-parabolic shape, but the power law with exponen¬ 
tial cutoff hypothesis is slightly preferred. The cutoff 
energy is (22 ± 10 s tat) GeV, while the spectral index 
below the cutoff is 1.46 ± 0.18 s tat* The source is not 
significantly detected at low energies. Using the same 
spectral binning as for the Cygnus Loop, the flux is 
measurable (TS^ n > 10) between 1 GeV and 72 GeV, 
making this source a candidate very-high-energy emit¬ 
ter. We also performed an unbinned likelihood analy¬ 
sis in time intervals of 60 days. The source is detected 
with TS 2 month > 10 in 13 out of 37 such intervals. 
This hint of variability supports the association with 
the background AGN, 2E 2051.7+2911. 


4. CONCLUSIONS 


Due to its proximity and angular size, the Cygnus 
Loop permits spatially resolved studies of different 
parts of the SNR, that interact with different com¬ 
ponents of surrounding medium. Particularly, it is 
known that the NE of the shell interacts with rela¬ 
tively dense medium and is thus bright in X-rays and 
optical emission lines compared to other parts of the 
shell. These inhomogeneities are likely to happen in 
other remnants from core-collapse supernovae, while 
remaining unnoticed due to lack of resolution of the 
instruments. Understanding the physical mechanisms 
that power the gamma-ray emission of the Cygnus 
Loop, and the differences between different regions of 
the shell, may help understand the variety of spectral 
shapes that SNRs display at gamma-ray energies. 

The fact that most of the gamma-ray emission from 
the Cygnus Loop follows closely the thermal X-ray 
emission from shocked matter supports the idea that 
most of its gamma-ray emission is emitted by interac¬ 
tions of accelerated hadrons with the dense medium. 
In this case, a low-energy break is expected in the 
spectr um due to the producti on threshold of neutral 
pions [Ackermann et aL||2013 .This is the case for the 
energy spectrum measured in this analysis, which has 
a maximum around 0.6 GeV. However, the SW por¬ 
tion of shell (which is brighter in radio and fainter in 
X-rays), has a different gamma-ray spectrum without 
indication of a spectral break. The explanation for 
the different gamma-ray properties, and the related 
radiative processes of the two regions, including the 
two-SNR scenario, is under investigation. 
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